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Abstract 

O 

Radiative transfer theory is developed for the 1304A triplet of 

atomic oxygen which includes the effects of pure absorption by molecular 

oxygen and variability of the temperature governing the relative populations 

3 

of the ground state P_. levels. Starting from an essentially complete 

intensity-height profile for the dayglow, the theory is used to determine 

the initial source function or the production rate. From the solution, we 

are able to discuss quantitatively the important dayglow excitation sources. 

The solutions are shown to be somewhat sensitive to the following parameters: 

(1) the oscillator strength f, (2) the temperature governing the Doppler 

3 

line width, (3) the temperature governing the Pj levels and (4) 
absorption below about 130 km. The assumption of complete frequency 
redistribution is applied and the basis of its validity is analysed with 
regard to the given problem. 

A source is required in the Fl region whose strength is about 200 

3 

excitations per cm per sec at 170 km. The altitude profile for the required 
source agrees well with those calculated for excitation of oxygen by supra- 
thermal photoelectrons, 

tThis work is based on a portion of the thesis submitted by D. J. Strickland 
to satisfy requirements for the Ph.D. in Physics at the University of 
Pittsburgh. 

*Present address: Department of Physics, University of Florida, Gainesville, 

Florida. 
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1.0 INTRODUCTION 

This is the first of two papers which will discuss transport of 
the resonance radiation of atomic oxygen in the Earth's upper atmosphere. 

In the present paper we treat the dayglow and in the second we shall con- 
sider the aurora. The radiation to be studied is a triplet at wavelengths 

O O 

of 1302-04-06A, abbreviated hereafter as 1304A. To date, a number of 
successful dayglow experiments have been performed in which altitude profiles 

O 

of 1304A radiation from about 100 km to 600 km have been recorded on sounding 
rockets (Fastie et al. , 1964a; Fastie and Crosswhite, 1964b; Kaplan et al. , 
1964; Fastie, 1968; Heath, 1968). Several analyses of the early dayglow 
observations have been published (Donahue and Fastie, 1964; Tohmatsu, 1964; 
Donahue, 1965; Tohmatsu, 1965) indicating that the observed height profiles 
and apparent emission rates require the presence of a source in the thermo- 
sphere more important for excitation of the triplet than absorption of the 
solar emission lines. These analyses, however, were based on observations 
carried only to Aerobee altitudes, i.e. , about 200 km. Such data do not 
extend high enough to show the form of the variation with height above the 
altitude of maximum brightness. This limitation, coupled with uncertainties 
about the absolute brightness of the alrglow lines and about the absolute 
flux and spectral profile of the solar lines, has left a certain ambiguity 
in the assessment of the relative importance of the internal and external 
sources of excitation (Donahue , 1965) . At the time these analyses were 
made, some high altitude data were available (Kaplan et al. , 1964) but are 
difficult to explain in terms of reasonable models and are suspected of 
background contamination. 

Recently, Fastie(1968) has observed the altitude profile of this 
dayglow feature from about 220 km to about 600 km during the flight of a 
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Javelin sounding rocket from Wallops Island in November 1964. The profile 
obtained at that time can be extended downward by normalizing it to the 
earlier Aerobee observations between 100 km and 220 km. Thus completed, 
it gives a sufficiently complete representation of the "typical" spatial 
profile of the airglow radiation field to reveal the distribution of the 
sources required to create it. This paper will show how a radiative trans- 
port calculation can determine these sources from the airglow information 
as well as how it can put some interesting limits on the range of atomic 
oxygen model atmospheres compatible with the observed radiation profile. 

The results require the presence of a strong internal source centered near 
170 km to excite the 3s^S^ level of atomic oxygen at a rate of about 200 
cm~3 sec - *. Regarding the density models, the variation of the resonance 
line brightness with altitude above 400 km where sunlight should be the 
principal source of excitation appears to demand more atomic oxygen in the 
atmosphere than the models appropriate for the time of the Javelin flight 
would call for. These points will be discussed in considerably greater detail 
after we have developed the transport theory for the 1304 triplet. 

2.0 Atomic Properties and Temperature Dependence 

The 1304 triplet arises from the transition (2p)^ ^ q ” (2p.)^ 

1 3 

(3s) which is designated as a resonance transition both because it is 

allowed and because the ground state (%) is involved. Such transitions 
are characterized by large oscillator strengths (f values) and in turn by 
large cross sections. The f value we originally chose for this work has 
the magnitude .021 (Parks et al. , 1967). With such a small oscillator 
strength it proved impossible to reconcile the experimental altitude profile 
with reasonable oxygen abundances. A much larger f value seems to be needed. 
Recent measurements by Lawrence 0-9^9) yield a value of 0.46 and F. Kaufman 
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(personal communication) reports that the value reported by Parks et al, 

(1967) is too low by a factor of about two. Theoretical f numbers have 

been calculated by Garstang (1961) (f - .051) and Kelly (1961) (f = .025). 

Using Lawrence’s ( 1969 ) value of .0 k6, we have recalculated our solutions and 

favor these in giving a physical interpretation for our results. As a point 

of interest, we present solutions in Section 6 for both f values showing the 

effects of changing this parameter. 

The f value is related to the line center cross section a for 

o 

absorption of resonance radiation by the relationship 

ire^ f z, s 

a = — v "• ■ > (l) 

o me AVp/m 

where the line has a Doppler profile and Av D is the Doppler width given by 



We see by Equation (2) that is a function of the kinetic temperature T. 

This temperature varies from about 200° to as high as perhaps 1000° over 

the altitude range of interest in this study. Selected values of a as a 

o 

function of kinetic temperature for f = .021 and .0 b6 are shown in Table 
1. The values shown have been used to obtain the results which will be 
presented in Section 6. 

The frequency dependent cross section a(x) is given by 

a(x) = o q <J>(x) 

where $ (x) is the line profile and x is a measure of frequency in units of 
Doppler widths from the line center 
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x » (v *- v )/Av 
O D 

In past work (Donahue and Fastie, 1964; Tohmatsu, 1964; Donahue, 1965; 
Tohmatsu, 1965) a pure Doppler line shape has been assumed. In this case 

2 

<Kx) * © X (3) 


It is a well known fact that the contributions from the natural wings of 
the line become important in case considerable multiple scattering or im- 
prisonment of the photons occurs. Such is the case for us in the lower 
regions of the atmosphere. Consequently, we have generalized <j>(x) to the 
form 


<f>(a, x) 


i f *~ y 

1T _ i a 2 + (x - y) 2 


dy 


(4) 


In this form of ij>(x) (the iun-normalizedv VMgt function) a is a parameter 
specifying the strength of the wings of the line, 

a - y . ■ 

4ir AVp 

Y being the transition probability or inverse lifetime of the excited state. 
The parameter a, here referred to as the natural damping coefficient, may 
also be related to f by the relationship 

a _ 2 8 1 Tre 2 f (5) 

X 2 g 2 me Av d 
o 

(see Mitchell and Zemansky, 1961; Equations 33 - 39) where g^ and g 2 are 
the statistical weights of the lower and upper states (g^/g 2 = 3 for the 
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oxygen triplet). In addition to the cross section a Q , selected values of 
£ as a function of temperature for f =* .021 and .046 are given in Table 1. 

It can be seen that <j>(a»x) depends on the kinetic temperature 
through both x and a. In a situation in which Av^ varies by more than a 
factor of two, as it does for the temperature range we encounter, it would 
be desirable to include the variation in the transport problem. To do so, 
however, would make the calculations considerably more difficult and costly 
to perform. As is generally done in a problem of this nature, we have used 
constant values throughout the medium but have varied T as a parameter in 
such a way as to span the expected temperature range. From the results we 
can estimate how much effect the isothermal approximation has on our solutions. 

The relative populations of atoms in the three levels of the 3 P 
ground term are assumed to be given by 

-AE /kT/ -AE./kT 

' g j e / l V <6> 

where j ■ 2,1,0 and where the statistical weights gj are 5,3 and 1 for the 
respective levels 3 P 2 » ^1 ^ p 0° We see that Pj is a function of tem- 

perature as are o 0 and 4(&,x). Table 2 gives p^ for various temperatures 
including the limiting values p^ * 5/9, 3/9 and 1/9 as T + », It should 
be pointed out that this temperature (hereafter referred to as the 
Boltzmann temperature T_) does not necessarily equal the kinetic temperature 

D 

governing the Doppler width, i.e., Tp. If the relative populations in the 

fine structure levels are influenced by collisions between atomic oxygen 

and ambient electrons (as suggested by Dalgarno and Degges, 1968), then Tg 

will differ from T^ and in fact will be greater than T^. Transport processes 

involving both the 1304 triplet and the far infrared radiations resulting 

3 

from the transition between the ^2 1 0 ^ eve ^- s can also influence the 



6 


relative populations. Their effect would be to increase the populations 
3 3 

of the and Pq levels and the distribution would then not be describable 
strictly in terms of a temperature. Because of collisions between atomic 
oxygen atoms and surrounding particles throughout most of the medium, we 
shall assume that transport processes have a negligible effect on the 
relative populations and that the populations can be represented by Tg which 
may, however, be greater than Tg. It is a much simplier task to treat the 
case in which Tg varies through the medium than the case in which Tg 

varies. Accordingly, we have considered the effect of a variable T„ in our 
calculations. Usually T B has been permitted to assume the ambient 
temperature but we have also considered an isothermal case with Tg ■ 1000°, 
chosen to represent a situation in which there is appreciable cooling of 
electrons by atomic oxygen, 

th 

With p^ a function of temperature, the optical depth for the j 
line at altitude z is given by 

00 

T j “ a o ^ Pj( z ') nCz 1 ) dz* (7) 

z 

As a function of frequency, 

Tj(x) » <j»(a,x) 

When the medium is isothermal, the depth for one line can be related to that 
of another by 


T 


k 


P J 
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Such a relationship permits the transfer equations to be reduced to forms 
having explicit dependence on just one of the three depths and effectively 
reduces the treatment of the transport to a two level problem (Donahue, 

1965) . For the non-isothermal problem, a different variation of p^ for 
each j forces us to deal explicitly with each of the three depths in the 
transfer equations. Computationally, however , this introduces no difficulties 
provided we accurately relate one depth to another at each point in the 
medium. 

© 

Because molecular oxygen begins to absorb the 1304A triplet below 

120 km in a significant fashion, this absorption process has been 

incorporated info the transfer equations. In fact solution of these 

equations including true absorption has been obtained only for the case of 

auroral excitation where it is of particular interest (Donahue and Strickland, 

1969). From those results, however, we can accurately estimate the effect 

on the dayglow (Section 7) . The non-resonant cross section for is 

© 

nearly independent of frequency in the neighborhood of the 1304A triplet 
and is taken to be 5 x 10 - ^ cm^ (Lee, 1955)* We designate 

the corresponding optical depth by 

f? 

t(z) = 5 x 10 j nCt^Xdz 1 
3.0 Transfer Equations 


3 . 1 General Approach 

The problem we set ourselves to is to determine the rate S (r) 

O A 
3 

at which atomic oxygen must be excited into the S^ level in order to 
produce the observed photon flux. The most likely contributions to S q are 
from electron impact excitation, cascading from higher states, and absorption 
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of resonance radiation of solar origin. Obtaining the solution is 
complicated by the fact that the oxygen medium is optically very thick. 

This means that a photon emitted by an atom excited by one of the primary 
excitation events will be absorbed and thus excite another atom before it 
has travelled very far from its original source. The same in turn will be 
true of the resonance photon subsequently emitted by this atom. The 
result of this multiple scattering or imprisonment of resonance radiation 
is that the steady state density of excited atoms is much larger than the 
value it would have if each photon emitted were to escape the oxygen medium 
without recapture. The unidirectional flux of photons at a given point in 
the medium is determined by this steady state excitation rate or source 
function S(jj). The rate of primary excitation S o (r) represents only a 
very minor part of S(r)» most of which is produced predominantly by capture 
of imprisoned photons. 

2 

The "apparent" column emission rate in photons/cm -sec-4ir ster 
observed at r in a direction cos \ from the vertical will be given by 


4irl(r,p) 

1 


S(r') T(r - r') dz'/y 

t.V- p-' Vi* 


( 8 ) 


where T(r - r 9 ) is a transmission function which depends on the number of 

* Of 

atoms between r and r* in the direction of observation. If we are given 

S-V 

4irl, the first step in determining the primary excitation rate is to 
invert Equation (8) and determine the final source function S(jf). 
Unfortunately the process can be carried out only if a model for the atomic 
distribution is first adopted in order to specify T(r - r 9 ). 

With S(r) and the model atmosphere specified, the initial source 
function S Q (r) can then be found since S(r) must satisfy an integral 
equation of the type 
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s <£> “ s 0 ( ^ + 


s <£> G <fc- r')dr' 


(9) 


where G is an appropriate kernel relating the emission rate S(r\) at r' 

MS T& 

to the excitation rate at r. The integral goes over the entire medium. 

* 

To solve Equation (9) is not a trivial exercise, because S q is usually very 
small compared to S and because the kernel is singular for vanishing values 
of its argument. As previously discussed, other complications also arise 
in practice because of the presence of pure absorption of the photons, 
notably by Oj, because variations in temperature cause variations in the 
width of the resonance absorption line profile and in the population of the 
fine structure levels of the ground state of atomic oxygen and finally 
because complex changes in the frequency of the photon occur upon scattering. 
All of these effects will be considered (not always rigorously) in the 
present treatment. 


3.2 Complete Frequency Redistribution 

In the initial formulation of the transport problem, a two level 
atom will be assumed. Accomodation of the fine structure splitting of the 
ground state of atomic oxygen will be made in the next section. It will 
further be assumed that the geometry of the medium is plane parallel with 
the density of absorbers functions of z only. This is a very good approx- 
imation in the present application where the sun is reasonably high in the 
sky and the scale height of oxygen is of the order of one percent of the 
radius of the Earth. 

In this section we shall investigate the criteria which determine 
the applicability of complete frequency redistribution (CFR) to a trans- 
port problem. By definition CFR requires that the frequency profiles for 
both initial excitation and scattering from a unit volume have the same 
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frequency dependence as the absorption profile. In any state of equilibrium 
in which the collision rate between excited atoms and surrounding particles 
exceeds the transition probability for the given state, the above condition 
is at least approximately satisfied; in particular, both the emission and 
absorption profiles are represented by the Voigt function. When the 
collision rate is significantly less than the transition probability, a 
unique relationship exists between the frequency and direction of an incoming 
photon and those of the outgoing photon for a given absorber of velocity 
v. in such a case if is not nearly as simple to establish that CFR is 
applicable. We shall now consider in a semi-quantitative way, how accurately 
the emission profiles can be represented by Voigt functions. 

For a plane parallel geometry, the basic equation of radiative 
transfer is 


PjT % - k(x,z) I(x,p,z) + e (x,p,z) (10) 

2 

where I is the frequency dependent intensity given in photons/cm -sec- 
ster-Av D in the direction p = cos 8, k(x,z) is the absorption coefficient 

3 

and e(x,p,z) is the emissivity or number of photons emitted/ cm -sec-ster- 
Av D . The absorption coefficient k(x,z) when expressed explicitly in terms 
of its frequency is k (z)$(a,x) , 4>(a,x) being the Voigt profile in this 
work. Expressed in terms of the optical depth t , Equation (10) becomes 

p4~- * - I (x , p , x) + S(x,p,x) (11) 

x 

where dt x « k Q <{>(a,x)dz and % S(x,p,x) , the total source function, is the ratio 
of £ to k(x,z). The intensity I is the same function whether we index it 
by t or z. For CFR, we note that S is not a function of x because e possesses 
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the same frequency dependence as does k. The formal solution to Equation 
(11) is 


Kx,y,z) - I (x,ii)e" T * (a » x) ' + S(x,M,T')e"l T ~ T 'l <t,(a ’ x)/w d[T'<|,(a,x)]/y (12) 


The first term on the right contains the contribution from the attenuated 
flux of radiation incident on the medium from outside. The second term 
represents the flux resulting from emission with the medium along direction 
y. A complication in using this solution arises from the fact that e is a 
function of I x not only at frequency x hut at all other frequencies as well. 
For the condition of radiative equilibrium, we may write e as 


e(x,8,z) ■ c Q (x,n,z) + 


l(x' ,8' ,z)R(x,x' ,8,8’ ,z) dx'd8' 


(13) 


This relationship expresses the condition that the emissivity must equal 
the rate at which photons of frequency x are created in a unit volume plus 
a scattered contribution from that volume . The quantity R, referred to as 
the redistribution function describes this scattering and is the probability 
for a unit volume that an incoming photon of frequency x' and direction 
ft ' is absorbed and re-emitted at frequency x in direction Q . For the situation 
in which collisions are unimportant , re-emission is isotropic, and the wings 
of the spectral lines are included, R has the form 


k (z) 


(x - X* ) 


sin0 


CSC 



sec 


0 . 

2 



ilk) 


exp - 



12 


where 0 is the angle of scattering. This form of R was originally introduced 
by Henyey (1940) and more recently analyzed by Hummer (1962) . A property 
of R which will be useful below is described by the following requirement 
of reversibility 


I 


R dx* 


k 

47 Ka , x) 


(15) 


If we wish, we can now immediately integrate Equation (15) over direction 
and obtain 


^ R dx* dO’ s» k o <Ka,x) 


(16) 


Because of the symmetry, we can also, if we wish, consider x' and 9 ' as 
either the incoming or the outgoing variables , If they are interpreted as outgoing 
variables, we must, in fact, obtain k Q ^ in Equation (16) since from the 
above definition of R, an integration over all outgoing directions and 
frequencies must leave us with the probability (per unit length) of absorption 
at frequency x . 

For the case in which collisions are important, the redistribution 
function, defined similarily to R above and here designated by R c has the 
simple form 


k 

R c (x,x') =^^f(a,x) <Ka,x') 


We immediately see that R c satisfies the same relationships as those expressed 
in Equations (15) and (16) for R and has the same symmetry properties. 
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Let us now consider Equation (13) , i.e. radiative equilibrium, when 

the intensity is independent of frequency although not necessarily isotropic. 

Equations (15) and (l6) show that the scattering contribution to e is just 

J(z) k <J> (a,x) where J(z) is the effective flux entering the unit volume 
2 

in photons/cm - sec - Because R c satisfies Equations (15) and (l 6), 

this result holds for R c as well as R even though each redistribution 
function describes entirely different scattering processes . A close study 
of the properties of the redistribution function R shows that when x' is small 
the frequency x of the re-emitted photon has a large probability of being 
anywhere in the frequency range where the Voigt profile is almost Doppler 
in shape - in the so called core of the line. When x* lies in the range 
of frequencies so large that the Voigt profile almost has the natural line 
shape the scattering is almost coherent and x is very nearly the same as x’ . 
Physically this behavior reflects a contest for photons between the .natural 
line absorption profile and the Maxwell-Bolt zmann distribution. For 
frequencies close enough to the line center the preference is for capture 
at the center of the natural line by sufficiently fast atoms while for frequencies 
far from the center capture in the natural wings by slow atoms predominates 
over capture by atoms in the tail of the Maxwell distribution. Thus, if 
one actually applies R to a transport problem, frequency redistribution 
occurs at small frequencies but gives way to coherency in the wings. If 
we now consider R c , we see that a photon absorbed at x' = 0 has a distinct 
probability - varying with the line profile <)>(x) - of being re-emitted in 
the wings, even the far wings. On the other hand, a photon captured far 
in the natural wings has its greatest probability of re-emission at x = 0. 

In the problem we treat in this paper the redistribution of frequencies 
is really described by R(x,x' ,r,r' ,z) . It is our hope, however, that throughout 
most of the medium, the intensity varies weakly enough over the important 
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frequency range that scattered radiation is roughly described by a Voigt 
profile, thus allowing us to replace R by the much simpler function R c . 

We emphasize, however, that R c is not the proper redistribution function 
physically, in our case. 

With complete frequency redistribution, the scattered radiation from 

a unit volume possesses a Voigt profile. From Equation ('13) » it then follows 

that the total emission rate e also possesses this profile provided the 

initial rate e does. To apply CFR validly 
o 

1) e Q must approximately possess a Voigt profile 

2) l(x,r,z) must be approximately independent of frequency. 

Regarding condition (2), it will usually be sufficient if the frequency 
independence extends out to about X s 3. Were l(x,r,z) independent of frequency 
over an unlimited range, we would find that most of the absorption occurs 

over the region |x| ^ 3 for the values of the damping coefficient appropriate 
to this study. 

In the case to be treated in this paper the initial source function 
should possess a profile closely resembling a Voigt profile because of the 
nature of its excitation. The dominant source of excitation is probably 
electron impact on atomic oxygen and thus will produce an initial source 
with roughly the required profile. To investigate the second criterion, 

O 

we have calculated the intensities of the 1302, 04, 06A lines at various 
depths in the medium as a function of both frequency and direction. 

Figure (l) shows examples of these at a depth = 600 for one of our 
earlier models with f = .02. The source functions from which these 
profiles were derived were themselves obtained in the CFR approximation. 

If the intensity implied by these source functions were to violate the 

second criterion seriously that in itself would have shown a gross inconsistency 

in the application of the approximation. We believe, however, that these results 
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present a realistic picture of the true intensities. This is based in 
part on several more critical and detailed analyses of the condition of CFR 
(see e.g. Holstein, 1947, Hearn, 1964, and Doschek, 1969). Holstein (1947), 
for example, was the first to show that interior to the medium, CFR will 
generally be valid in the core of the line. 

Considering Figure (1) , we see that the intensity profile is 

O 

generally flat to about x * 2, The reason the 1306A profiles are the 

narrowest is that the total effective optical depth t q is smallest for the 
© 

1306A line. Although the frequency independent portion of none of these 

profiles extends to values of x as large as we would like, say to |x| > 3, 

still an appreciable fraction of the absorption which would result from a 

perfectly flat profile does come for values of x between +2 and -2. For 

the intensity profiles shown, the resulting emissivities for scattered 

radiation would closely resemble the Voigt profile in the core (or the 
2 

e portion of the line) but would show a deficiency in the wings. And, 
indeed, in the real case described by the exact redistribution function, 
the emissivity is deficient in the wings, because of the tendency of photons 
in the wings to escape the medium after a few scatterings. 

The deficiency of wing photons illustrated in Figure (1) shows, 
however, that the CFR approximation unrealistically enhances the emissivity 
in the wings of the line. If we now consider Figure (2) which shows 
intensities evaluated very close to the upper boundary, we see that CFR tends 
to break down there. It is not very far from the boundary, however, that 
the situation regarding the intensity profiles improves considerably. So, 
over most regions of the medium, as long as the examples in Figure (1) 
are typical of the true intensities, CFR should not yield results 
significantly different than those obtained in an exact calculation. The 
investigations referred to above all indicate that this is indeed the case. 
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Now assuming CFR, the total emission rate can be expressed as 
e(x,ft, z) =« e(z) <j>(a»x) /4ir«^T 

where <j>//ir is the normalized Voigt profile and l/4ir normalizes e for 
isotropic emission. The source function, originally introduced as S(x,y,x) 
becomes 


\ e(z) 1 
s <x,p,T) «— 477JT 

o 

and is seen to be independent of x and y, In the treatment to follow, we 
choose to define the total source function by 


SCT) - 

O 

2 

which represents the total emission rate per cm and per unit optical depth. 
The total emissivity is then 


e(x,z) 


k 


o 


S(t) 


4ir/ir 


(17) 


Similarily 


e Q (x,z) 


k _ S (t) 
o o 


<Ka,x) 

4ir/ir 


(18) 


Replacing, R by R c in Equation (13) and dividing out the common factor 
k i/UWit we now obtain 

S(x) - S q (t) + 


I(x’ ,y' ,xH(a,x’)dx’ 2Trdy' 


(19) 
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Equation (12) allows us to express this as 


S(t) - S (t) + 
o 



T~ T 


auj. 

M' 


( 20 ) 


Here, the contribution of the attenuated external radiation source has been 

lumped in the S q term along with internal excitation sources. It should 

be pointed out that if one wishes to solve this equation by initially 

specifying S q for an external radiation source, special considerations must 

be given as to how S is evaluated. Deep in the medium if S were to 

o o 

represent the actual rate of initial excitation produced by the external 

radiation source we would find it heavily dominated by wing photons due to 

the ever increasing reversal of the attenuated external flux as it penetrates 

the medium. Most of the photons resulting from this excitation are far 

enough in the wings that they would escape the medium almost immediately. 

Hence, to choose an effective CFR S , we would have to judiciously discard 

o 

most of actual wing dominated S q , discarding an ever increasing fraction 

deeper into the medium (Donahue, 1965). Fortunately, we can ignore this 

difficulty since our objective is to derive S from the intensity. 

o 

If we now introduce the transport kernel 
H( |t - t* | ) = 2 ^- <t> 2 el T " T dx' = ^ <() 2 E 1 ( |t- t‘ |jdx* 

J 

( 21 ) 

Equation (20) takes the form 

» 

S(x) * S q (t) + S(V) H(|t - t 1 | )dx 1 (22) 

4 

H is the probability per unit optical depth at x that a photon upon emission 
from the level x* will reach the level x and be absorbed. 
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The intensity resulting from emission within the medium can be 
expressed as 


I(x,y,T) S(t') • ^ a ^ ) 


-|t - t' |<j>/p 


a,x)dT , /y (23) 


Integration over x permits Equation (23) to be written in the form 


4itI(t s ji) * \ S(t') T(|x - x'l/p) dx'/p 


(24) 


where 


, . 1 0 -|t - t’ I <J>/p 

T(|t - x' l/u) \ $ e dx (25) 

The transmission function, T(jt - x'|/u) is the probability that a photon* 
will be emitted and will travel the distance |t - x' |/y without an absorption* 
Equations (22) and (25) are the CFR forms of the transfer equations for a 
two level atom where the medium can be represented by a plane parallel 
geometry* In the following section we shall modify this formulation by 
incorporating in it absorption and the fine structure of the ground state 
of atomic oxygen. 

3.3 Introduction of pure absorption and fine structure 

Pure absorption has the effect of increasing the effective 
attenuation along a photon's path. We take account of this loss of photons 
by properly altering the probability functions T and H. The new transport 
functions ,T'(|x-x'|/y, |t-t'|/y) and H' ( | x - t ' | , 1 1 — t * | ) become 
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t'(|t - t’|/u. It - t * | /m) 


i 

TtT 


<j>(a,x) 



t' 1 4>(a,x) + 1 1 - t* | ]/y 

dx 


It - t'| 
e T(|t - t’ |/ y) 


(26) 


and 

H'(|t - T 1 1, |t - t' |) - Yfe |^» 2 (a,x) E^Cl't - t' |<Ka,x) + |t - t' |)dx (27) 
To take into account the fine structure of the ground state of 

3 

oxygen we note that all emissions arise from the single level S^. This 
allows us to designate S as the total source function despite the fact that 
the resulting emissions separate into three distinct lines. If, however, 

O 

we desire that part of S which produces the 1302A line, for example, we 
have 


S 


2 



S 


or more generally for the 



line 


« gj ’ S (28) 

where g ' is the ratio g./Z g, . The same considerations also apply to the 
j 3 k 

initial source function S . In solving the transfer equations, S and S 

o o 

may be expressed as a function of any of the three depths. He choose to 

write them as functions of T 2 and as such they are then given in units of 
2 

photons /cm -sec-unit T 2 » If we wish to transform from one depth to another, 
we can use the transformation 

S< V - J S( V 


(29) 
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and do similarity for S Q . The contribution of S to the j** 1 line (Equation 
28) may then be written as a function of any chosen depth x^ in the form 


yv - y s( y 


(30) 


til 

where k need not equal j and (x^) is the rate at which the j line is 
emitted per unit x^- 

Now, considering the first of the transfer equations, the intensity 
in the j** 1 line is 


4irI^(Xj,p) 


yy )r 




|t * t' |/p)dx^'/y 


(31) 


Experimentally, only the sum of the three lines was measured in a 
vertically upward direction and as such the observed intensity is expressed 
by 


4irl(x 2 ) 


E g, \ S(x j ’)TV(|x j 


Tj'l, It -»t , |)dT j 


(32) 


where x 2 on the left is used simply as an index of the level at which the 
intensity is evaluated.. 

The total source function S(x 2 ) is comprised of S q (x 2 ) plus 
scattered contributions arising from excitation by each of the three lines. 
For the j th line, this latter contribution is 


S(t j , )B f (.| 


- T , 


j J 


t - t' | )dXj 


This term, however, is in the form of an excitation rate per unit depth x^ 
in the appropriate line. To add the contributions of the three lines, they 
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must each be changed to rates per unit x 2 and from Equation (29) , this 
requires that each be multiplied by p^ /p 2 . The transport equation then 
becomes 


s(* 2 ) 


w 


+ E 

3 


if. 


£i 

3 P 2 


S<T j )H(|T j 


- T , 


ft - t' 


l>*V 


(33) 


From the nature of our problem, considering the large optical depths involved 
plus the addition of pure absorption, difficulty was encountered in obtaining 
a smooth function for S q . Because solving the above equation required much 
of our effort, a discussion of the sources of error and the final workable 
form for this equation as well as Equation (32) are given in the Appendix. 


4.0 Data and Atmospheric Models 
The -height profile of the intensity is shown in Figure 3 and 
represents a composite from two experiments (Fastie, 1968). During these 
flights the solar zenith angle was <|0$. For the time of the high altitude 
flight, the Jacchia (0965$) model predicts an exospheric temperature of 850°. 
An atmospheric model appropriate to this temperature is given in Table 3 
(A. I. Stewart, personal communication). As will be discussed shortly, a 
requirement regarding the behavior of the total source function S(x) suggests 
the possibility that a higher exospheric temperature existed during the time 
of the high altitude flight. To account for this possibility, we have also 
considered a model having ancexospheric temperature of 1000° (Table 4). 
that part of the model above 120 km was taken from U.S. Standard Atmosphere 
Supplements, 1966 and that below from Colegrove et al. (1966), the latter 
being scaled so as to join smoothly at 120 km. In Tables 3 and 4, the 
optical depths x 2 , tj and " T q are included to show their range in the 
problem. The particular sets given refer to T D « 85b* and 1000° respectively 
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and Tg equal to the appropriate ambient temperature profiles . 

5.0 Total Source Function Near the Upper Boundary 

In the course of solving the transfer problem, we have obtained 
solutions for both f = .021 and f * .046. The basic model used with f * .021 
has a 1000° exospheric temperature (Table 4). Two models were used with 
f » .046, the 1000 model just mentioned and another with an 850 temperature 
(Table 3). Table 5 shows the particular models for the various solutions 
we have obtained. 

Based on the Jacchia model, we originally solved the transfer 
problem using the 850° model and f = .021. Trouble immediately developed 
because of an upturning of S(x 2 ) nfear X 2 “ 0. A requirement inherent in 
the transport theory is that $(x) monotonically decrease as the boundary 
is approached and, in fact, possess an infinite slope at x * 0. The con- 
dition at the boundary is readily demonstrated if we differentiate the 

i 

transport equation. Initially we have 


dS 

dx 





)dx' 


(34) 


From equation (21) , 


•4H , and consequently , by replacing -g=- with 
dx ' 1 


Hx"* and integrating by parts, we obtain 


dS 

dx 


dS 

m + S(o)H(x) - S(x q )H(x 6 



(t’)H(|x 



|)dx' 


(35) 


As x •* 0, H(x)-*- 00 which is sufficient to establish that ». Away from 
the boundary S(x)will continue to increase to a peak located well into the 
medium. Even with an external source of excitation, S(x) peaks well away 
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from the boundary because of imprisonment which is strongest at the center 
of the medium. 

To correct the behavior of S(t 2) requires . greater optical depths 
near the 600 km level. For this reason, the 1000® model was applied to the 
problem and this resulted in a considerable improvement although a slight 
upturning still persisted. Another method of increasing the optical depth 
is to increase the f value. Theref Qre , we adopted Lawrence* s (1969) value 
of .046 and recalculated S(x 2 ) for both exospheric temperatures. The 
results are shown in Figure 4 along with two other solutions to be discussed 
shortly. The positive gradient in 5 still persists at 850° but is 
essentially or nearly corrected at 1000®. 

If we are to believe that the 850® model is actually the more appro- 
priate one, we must then consider sources of this pathological behavior 
in S other than a lack of optical depth. One such source might be the 
breakdown of CFR near the boundary. This can be ruled out as a problem 
however, since for small optical depths, the CFR solution S(t) must, in fact, 
approach the true S(x) and both in turn must approach the solution for the 
case of no attenuation. To represent this latter solution near the boundary, 
we divide the medium into zones and consider average values of S over each 
zone. For the first two zones, 

* 4itI 1 /At 1 


and 


S 2 :« (4itI 2 - Airlj) / At 2 

where 4irl^ and 4 ttI 2 are the observed intensities at the bottom of zones 1 
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and 2 and Ax^ and Ax£ are the zone widths <> The resulting solutions for the 
850® and 1000® models are included in Figure (4}> along with the CFR solutions . 
Because the actual emission profile corresponding to the true S(x) shows 
a predominance of wing photons near the boundary, initially the true solution 
will not differ much from that with no attenuation but more specifically 
will lie between the solutions for this and the CFR case., Since all solutions 
exhibit an upturning in Figure (4) , so must the exact solution and we thus 
conclude that the assumption of CFR is not responsible for the effect. 

Another possibility is that 4x1 actually decreased more rapidly 
with increasing height than the data indicated. This is difficult to 
believe, however, since the unresolved triplet was observed with a 
spectrophotometer, the resolution of Which was sufficient to rule out 
spectral contamination. Unless there was a time dependent source of 
absorption not properly corrected for, the relative intensities recorded 
should be reliable. Although we have not satisfactorily resolved the 
question, the best explanation seems to be that larger optical depths are 
called for. If we assume the f value is no larger than .046, this requires 
us to abandon the 850® model in favor of one with a higher exospheric 
temperature, one in the vicinity of 1000° , or somehow to introduce considerably 
more oxygen at high altitudes. 

6.0 Temperature Effects Exhibited by the Transfer Solutions 

In Table 5, we have listed models in terms of the parameters f, 

T_ and T_. The solutions we have obtained for them show sufficient 
D iJ 

differences that a brief analysis of the various effects is in order. Before 
proceeding to compare the various solutions, it is worth noting that a 
change in f has exactly the same effect as changing T^. frovided T Q is taken 
to be constant throughout the medium, then T^ as well as f enter the transfer 
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equations only through and a (Equations (1) and (5)). The relationship 
is an inverse one such that increasing f has the same effect as decreasing 
Tjj by an appropriate amount. In the following analysis, we first show the 
effects produced by changing T g and then T^. It will then be straightforward 
to compare solutions for differing f values since the effect is the same 
as changing Tg. 

To show the effects of changing T„ and T_, we consider solutions 

o U 

for f » .046 and T * 850®. Changing the exospheric temperature from 850® 
to 1000® has a negligible effect on at low altitude and so we shall 
confine this discussion to the 850® model. In figure (5) are shown the 


total source functions S 


3V (b ) 


and Sg ^ ese are solutions for 


a kinetic temperature of 300®K but for a variable Tg and for a Tg constant 
at 1000® respectively. The subscript b denotes the models in which T is 

6X 


850® and f is 0.046. Because each of these solutions relates to the same 


intensify profile, if follows that model 3C is the more opaque model 
requiring a greater source function than 3V. If we consider individual 
optical depths, we see from Equation (7) that and Tq of 3C increase 
relative to those of 3V as the altitude decreases while the reverse 
happens to The fact that Sg c is greater than indicates that relative 

increases in and Tq more than compensate for the relative decrease in 
t 2 ° Thus the effective opacity of 3C is greater than 3V as measured by 
the transmission function relating S(z) to 4irl. 

The initial source functions S (z) as obtained from S(z) by use 

O 

of Equations (A4) are compared for models 3C and 3V in Figure (6). Generally 

the medium with the greater opacity will provide greater amplification of 

S q by imprisonment; as the escape function E, is smaller. Thus a smaller S Q is 

required to produce a given S(z) at a given altitude. This effect is not 

3C 3V 

great enough, however, to reduce S q very much relative to S q despite 
the greater opacity of 3C in comparison with 3V. 
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The effect of changing Doppler temperatures can be seen in Figure 
(6) where S(z) are compared for models 3C(b) and 8c(b), in which the Doppler 
temperatures are 300° and 850° respectively but T_ is constant at 1000°K. 

15 

The line center optical depths are larger for the 300° case than for the 
850° case by a factor of 1.68, Thus where transport occurs mainly in the 
Doppler core of the lines model 8C is more transparent than model 3C and 
requires a smaller source function. This is the case above 160 km. At 
lower altitudes, however, the optical depths become so great that the 
wings dominate the transport. In this situation, the medium of 3C actually 
becomes more transparent than that of 8C thus requiring less source function 
as Figure (6) indicates. Quantitatively this effect can be understood if 
4>(x) is represented, as it can be, by 

j 1 ^ —x*~ 

fix) = e , x « x c 


♦ <*> - 


a 1 


x 


X*>> X 
2 C 


where 


-xf 
e b 


a 

77 


(36) 


Then 


+00 


T(|t - T * | ) « 


1 

/r 


2 - T 

—x e 1 


- t' 


t ' e 


-x 1 - 


dx 


(37) 
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when |t - t'| is small, and 


T(j 


*'l> 


s 2a 

~ IT 


1 - j T - T* 

2 * 


a//ir x 


d?c 


(38) 


when | t - t ’ | is large . Between the 120 km and 140 km levels the optical 

© 

distances in the 1304A line are 7800 and 4637 for the 300° and 850* Doppler 

temperatures respectively. However, the transport functions, given practically 
, -4 -4 

by Equation (38) are 8.39 x 10 and 8.00 x 10 for 3C and 8C respectively. 

The initial source functions S Q (z) required to produce S^q and 

3C 8C 

Sgg are shown in Figure (8). S q dominates S q relatively more than 

dominates Sg c even at high altitude. This indicates that medium 3C is 

not particularly efficient compared to 8C in trapping photons despite the 

greater line center opacity. The reason is that large slant paths are 

Involved in determining the escape probability E and the transport kernel 

H. Because of the effect of the wings the model 8C remains a more effective 

imprisoner of photons than 3C up to fairly high altitudes. 

In Figures (9) and (10), in addition to the initial source 

functions just discussed for models 3C(b) , 3V(b) and 8C(b) , are shown the 

excitation rates required for models 8V(b) , 3C(a) , 3V(a) , 10C(a) and 10V(a). 

The last four cases, shown in Figure (10), are the models in which the 

oscillator strenth assumed was 0.021 and the exospheric temperature was 

10Oft*K*JThe similarity of the profiles for cases 8(b) and 3(a) is a 

consequence of the effect we have mentioned above, that increasing the 

oscillator strength is equivalent to reducing the Doppler temperature. The 

development of a bimodal structure in S q as the Doppler width Increases 

and the line center opacity d£MM*h’es. Is Interesting in view of the possible 

3 

importance of dissociative photon excitation as a source of atoms at low 
altitudes. We shall discuss this point at greater length in the next section. 
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7.0 Excitation Mechanisms 

In Figure (11) the initial source functions or excitation rates 

3 

for the 0( Sj) term required to account for the rocket intensities are 
shown for two models. These are 8C(b) ,(850® Exospheric temperature and 
Doppler temperature, constant 1000° Boltzmann temperature and an oscillator 
strength of 0.046) along with model 10C(c) which is the same as model 8C(b) 
except that T w and are 1000®. The second model therefore contains more 
atomic oxygen than the first at high altitudes and has a considerably 
greater line center opacity. In both models there are clearly two principal 
source regions. One source maximizes above 400 km and is about two orders 
of magnitude weaker than the other which has a maximum at about 170 km. 

Presumably the high altitude source is provided by absorption 
of solar photons in the three solar 01 resonance emission lines. The 
contribution to S q from this source should be given by 

S o (z) - Z(irFv)^/w AVjj T(r k /p b ) o q p k n(z) 
k 

where ji q , the cosine of the solar zenith angle, is 0.5 in our case and 

(irFv)^ is the effective solar flux at the center of the k line in units 
2 

of photons/cm sec Hz. According to Tousey (1964) and Tousey et al. (1964), 
when AVp is appropriate to an 850® Doppler temperature the effective flux 
within Av D should be given by 

7 2 

(irFv^ /it AVp - 5.2 x 10 ph/cm sec 

(irFv)^ /ir Av D = 8.3 x 10 ph/cm sec 

7 2 

(ttFw)q /it Av^ ■ 10.4 x 10 ph/cm sec 
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These values assume that the solar lines are square in profile and 0.2A 

wide. They are thus lower limits unless the lines are seriously reversed. 

In case of reversal we note that S would increase with decreasing altitude 

o 

more rapidly than we calculate here. 

In Figure (11) S o from this solar source is plotted for the 850°K 

atmosphere. For the 1000° model the maximum occurs higher than for the 

850°K model but has about the same maximum strength. The agreement between 

the predicted source shape and strength and that calculated from the air- 

glow data would be fairly good were the solar flux about one third as great 

as we have taken It to be, the airglow three times as bright or there were some 

combination of reduced solar flux and Increased airglow emission rate. The 

solar S also seems to be shifted about one scale height above the S required 
o o 

from the airglow data for the 850®K model. Increasing the oxygen content 
of the atmosphere does not tend to shift one curve relative to the other, 
but a reversal in the solar line profile could bring the solar source 
downward. 

Furthermore, it is well to recall that the assumption of CFR is 
more than questionable in the last few optical depths of a very thick 
medium. The outgoing radiation flux is seriously reversed (Figure 2) , 

The deficiency in the real case means that CFR, assuming as it does a flux 
independent of frequency overestimates the outgoing flux available for 
scattering near the line center. Hence it tends to underestimate the 
magnitude of S q really required to produce a given scattered radiation 
field. This remark should serve as a caveat for attempts to deduce the 
density of an atmospheric constituent from the measurement of the brightness 
of scattered resonance radiation. Even in the optically thin boundary of 
an atmosphere which becomes optically thick below it is necessary to know 
the line profile and absolute flux of the radiation emerging from the 
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interior of the atmosphere as well as the external radiation to obtain the 
abundance of scatterers from the scattered emission rate. 

The inadequency of the solar flux source in accounting for the 
magnitude and height profile of the 1304 airglow brightness shows clearly 
in Figure (11) and confirms the original conclusion of Donahue and Fas tie (1964) 
that an internal excitation source in the Fl region is required. The 
profile of the source required is shown in Figure (12) for the model 8C(b). 

This model is chosen as the one most likely to represent realistically 
atmospheric conditions at the time of the Javelin experiment. Following 
suggestions by Dalgamo (1964) and Barth (1963) that energetic photo- 
electrons should provide the necessary source there have been a number of 

3 

calculations of the rate of excitation of the term by this mechanism. 
(Tohmatsu, 1964s Green and Barth, 1967$ Dalgamo et al», 1969 )• 

In Figure (12) we show the rates calculated by Prasad (1969) and 
A. I. Stewart (personal communication, Dalgamo et al. 1969) for conditions 
similar to those assumed in our calculation. The Stewart calculation is, 
in fact, for the same model. The Prasad calculation is based on our 1000 a 
model, but allows for escape of photo-electrons while Stewart assumes local 
energy loss. Allowance for escape should bring the Stewart curve into 
fairly close agreement with the required S q above 230 km. The excitation 
rates have all been normalized at 170 km. That of Stewart in particular 
is more sharply peaked and reaches a maximum somewhat lower than the required 
source. However, in general the agreement between the calculated excitation 
rate by photo-electrons and the rate needed to explain the day glow data is 
good. At low altitude absorption by 0 o requires that our values for S be 
increased below 150 km (Donahue and Strickland, 1969). The correction is 
about 20% at 130 km and should remove the discrepancy to be seen in Figure 
(12) at low altitude. 
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On the other hand there are other potential sources of excitation 
which might be significant. Two of these are dissociative processes involving 
0^, one upon electron impact the other following photon absorption. The 

O 

threshold for the photo-dissociative source is at 840A and there is an edge 

in the absorption cross section at that wave length. The magnitude of the 

increase in the cross section suggests that the partial cross-section 

-17 2 • 

involving excitation might be about 2 x 10 cm at 800A (Cook and Metzger 

1964 ) • 

3 

The rate of production of the oxygen 3s term by this process 

has been calculated for both the 850°K and 1000®K atmospheric models 

treated in this paper. Solar fluxes and absorption cross sections have been 

taken from the tabulations of Hinteregger et al. (1965) and the cross section 

for photo-dissociation estimated from the measurements of Cook and Metzger ( 1964 ). 

-3 -1 

The results are Chapman layers with peak rates of 54 cm sec at 155 km 

—3 —1 

for the 1000® atmosphere and 70 cm sec at 145 km for the 850° atmosphere. 

e 

Only the flux down to 790A is effective in this process owing to the strong 
absorption of N 2 below 790A. A possible exception is the solar 304A line 
which penetrates to about 120 km. A portion of this line may produce 
photodissociative excitation of O 2 in an important amount despite a small 
cross section because of the large solar flux. However, the excited atoms 
should have large velocities and only a small part of the radiation emitted 
can be imprisoned. In fact the radiation emitted from a large fraction of 
the excited atoms produced by dissociating 0 2 may be Doppler shifted so much 
that it has little chance of being absorbed locally. Thus a CFR analysis 
of the intensity-height profile would not call for a local maximum in S q 
from such a source. Although the most transparent models such as 10(a) 
call for a bimodal source with a maximum near the level where the dissociative 
process is expected to be effective, such models seem to be sufficiently 
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unrealistic as to lead us (regretfully) to regard this interesting feature 
as fortuitous. 


8. Summary 

Applying the techniques of radiative transfer theory to the 

@ 

analysis of the 01 1304A dayglow, we have determined the initial source 

3 

function or production rate of the level of oxygen starting with an 
observed intensity height profile and a variety of model atmospheres. The 
calculations have taken into account the variation of the Boltzmann 
temperature determining the population of the ground state levels and pure 
absorption by 0^. Complete frequency redistribution and a constant Doppler 
temperature have been assumed although different temperatures have been 
considered. In both the case of Doppler and the Boltzmann temperatures, 
n on-negligible effects are observed as these are changed. 

By the nature of the calculations, the folution represents 
excitation by all mechanisms that were actually operating during the time 
of the experiment. The dominant sources of excitation appear to be photo- 
electron impact on atomic oxygen below 400 km and resonances absorption of 
the 1304 solar lines above this altitude, the former being, by far, the 
stronger source. Additional sources of secondary importance are cascading 

3 

into the S level from higher states of oxygen, radiative recombination and 
dissociative excitation of 0 ^. 
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Appendix 


To present the final workable forms of the transfer equations, 
we first consider the two level problem without pure absorption. Once 
the theory is thus formulated, we add pure absorption and finally take 
account of the fine structure of the oxygen ground state. Our numerical 
approach is to divide the medium into N variable zones, with the upper 
and lower boundaries of the n zone specified by the depths and 

+ y Variable zones are applied to take account of rapid variation of 
S(t) near the boundaries but otherwise slow variation over the rest of the 
medium. At the bottom of the n tn zone with optical depth the 

intensity has the form 


m + 1 


1 

4irl -IS 
n m 

m - n 


I<T n + 1 - T ' )dT ' 
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where S is the average value of S(t') over the m zone. In terms of the 
m 

so called width function 
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the intensity becomes 
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(A3) 


The matrix * ^ is triangular and thus inversion techniques are not required 

to solve for S . Beginning with the first zone we have 
m 
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Si = 4 jti 1 / $ 

followed by 

S 2 " (1WI 2 ” S j> 21^22 

and so on. 

Once S(t) is determined, the initial source function S q (t) can be 
found most readily if Equation (22) is written in the form (Donahue and 
Fastie, 1964) 

’ 

S q (t) * S(t) E (t) * J[S(t') - s(r)3 H( I t-t (A4) 

The function E(t), referred to as the escape function and defined by 

T o 

E(t) = 1- J H(| r-r» | Jd-r* (A5) 

o 

is the probability that a photon originating at t will escape the medium 

without being resonantly absorbed. The advantage of Equation (A4) is that 

S and the integral term have been modified such that their magnitudes are 

similar to S Q . In the form of Equation (22), an accurate solution of S Q 

for optical depths as large as those in this problem is not possible because 

of imprisonment which produces typical differences between S and S 0 of three 

or four orders of magnitude. 

In terms of N zones. Equation (A4) becomes 
T m+1 

S on" S A - “ f [S(T ' > - s n ] H( I V T ' I )dT ' (A6) 

m=l J 

T 

m 

where the functions S Qn> S n , E n , and H are referenced to t n , the center 
value of the n^ zone. We refer to the integral within the summation as 
the transport integral and designate it by 
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K = 


nm 


T m+1 

J [S(t’) - S n ] H(|f n -T'|)dT» 


(at) 


Most of the computational work involved in obtaining an accurate solution 
for S 0 is concerned with the proper numerical evaluation of this integral. 
Computationally, it is very useful to note that without pure absorption both 
E(t) and H(|t-t'| ) are related to the function 

f 1 

dp 


£ iV T 'i> = 2 


dx <J>(a,x)e 


- T n -T’ j <J>(a,x)/p 


(A8) 


which represents the probability that a photon emitted at t' will travel to 
T n without being absorbed. The relationships are 




( A9) 


and 


H dv T 'i) " - 


(A10) 


where * refers to t* > x . In terms of 6(1^ ” T '|) 


m+1 


K 


nm 


= - l [s( t * ) - si ax’ 


)dx' 


(All) 


m 


Integration by parts leads to the following two forms of K nm , the first for 
the case in which S(t') is assumed not to vary within the zone and the second 
in which it is assumed to vary linearly as a function of t ' . For S constant 



36 


L = [S m “ S J l^ (AT n J I 

nm m n n ,m+l n ,m 1 


(A12) 


while for a linear variation of S over each zone 


M “ t ^ AS n m+1 ^ ^ X r> m+1 ^ ” AS „ m £ ( AX nrn ^ 

nm — 1 n s m+x n ,m+l n ,321 nm 


(A13) 
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It is interesting to note that the term in large brackets in M depends on 

the properties of the model alone and not on the source function. 

Unless the zones are made so narrow and so numerous as to destroy 

the advantages for which the technique of zone division is designated, the 

approximation that S is constant for each zone can lead to serious error. 

Because H' decreases as x' varies across the zone, the effective average 

value of S(t') is found closer to tv than the center of the zone. The error 

n 

arises in the transport integral because is subtracted from S(t’) and the 
difference becomes very sensitive to S(t’) when S(t') - S q . Approximating 
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the real S(x') by a linear function over each zone considerably improves 
the accuracy of and this approximation has been used to obtain the 
results in this paper. In terms of the transport equation for S q may 
now be written as 


S 

on 


n 

S E + Z M 
n n , nm 


(A15) 


a form expressed in terms of the transport functions W(x) and £(x) . 

Two effects complicate these procedures for obtaining S(t) and 
S q (t) in the aurora and airglow. One is absorption of the resonance photons 
by molecular oxygen, the other is the triplet nature of the radiation. 

Keeping to the two level formulation of the theory for the present, absorption 
by Og can be taken into account when S(x) is deduced from 4irl(x) by letting 

a nm ba the °P tlcal d *P th £or absor P tl “ h ? °2 b,t ”“ n tha ce " tar the " th 
zone and x n + The transmission function in Equation (Al) then becomes 
T’(T n + ^ - t f , t^ + ^ - t f .) as defined in Equation (26). The result is 
that Equation (A3) is then changed to 


4irl 


I 

Z 


S e nm (W 


(A16) 


The effect is to increase the values of in regions where absorption is 
appreciable. 

Similarly, in the transport Equations (A4) and (A6) H’(|t - t'|, 

|t - t'|) as given by Equation (27) must replace H(|t - t'|). The transport 
integral (A7) becomes 

[S(t’) - SJ H'(|t q - T*|,|t n - t’|) dx» (A17) 

x 
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If t' is taken to be constant through each zone then 



T’1,1 Vt’D* 





t' |); 




(A18) 


where 


, ( | t - 
1 1 n 


x’, t 


nm 


1 

277 


,1 




-( j x 

dx $(a,x) e n 


L 


t’ j$(a,x) + t nm )/y 


(A19) 



4>(a,x) Eg[( |x n - t’| <f>(a,x) + t^)] dx 

<90 


Here t is the optical depth in 0„ from the center of strip n to the 
nm d 

center of strip m and E^ is the second exponential integral. 

If, however, we consider a linear variation in the t ' across a 
zone expressed by 


I? - t ' I * a. It - t * I + 8 
1 n ' t 1 n 1 t 


H' is defined as the divergence of where 


(A20) 


L 


;CF n - T' 


t n - f|, «,.) 


l 

277 


» 2 ( a >*> E (I? - T » 
+(a,x) + o t B 2 MT n t 


a,x) + |t n - t * | ) dx 
(A21) 


For m 4 n, the effect on K of replacing t by a variable absorption for 

nm am 

each zero is very small. Therefore, in practice, Oj absorption between zone 
n and m is taken into account by using the function cf ^ in evaluating the 
kernel H'. We need not even consider the case m ■ n since by the nature of 
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K , this term is virtually zero. Now L , the transport integral which 
nm nm 

applies to the case of constant S(t') in each zone, may be written in the 
form 




m + 1* 


t ) 
nm 


- £i< At 


nm’ 


Wl 


(A22) 


The second form M nm , appropriate in the case of a linear variation in S(t’) 
over each zone, becomes 







t ) 
nm 


m + 1 


4” 

a 


n 



(A23) 


after an integration of Equation (A17) by parts. After some manipulation, 
the integral term takes the form 


-t 


nm 


(w - w , ,) 

4 n,m n,m + 1 


(A24) 


“ t[ A t „ m + 1 £( AT n n, + 1» " At n 

z n,m +1 1 n,m t x nm n 9 w , x nm nm 


nm nm 


where 


€ 


^(t) 
nm nm 



Insertion of Equation (A24) in 

form for M 

nm 


, . <j>(a,x) + t ) - E 0 (At $(a,x) + t ) dx 
+ 1 r ’ nm 2 n,m * nil 

Equation (A23) leads to the final two level 



4o 


M 


nm 


± (AS . i f , (At ^ v , t ) - AS . (At ,t ) 
n,m + 1^-1 n,m + 1 nm n,m'— 1 n,m’ nm 

-t 

e nm 1 r 

+ “n'-T- «n,» ' M n,» + 1> * 2 «*„,» + 1 <£>*„,» + 1> W (A26 > 


- At (£,(At , t ))■■-£ (t )] 

n,m^-l n,nr nm ' nm nm 


which is to replace the earlier version in Equation (A13) . 

The escape function E must also be modified to accommodate pure 

n 

absorption. Beginning with the definition of Equation (A5) and using 
Equation (A10) this function becomes 

n 5 n+1 . r 

E n »l- I \ (±) r'|,|t n - t'|) dT' (A27) 

m ■ 1 J 


The sum is often very dose to unity and E n very small for an optically thick 
medium. The largest contribution to the series, even if At^ is of the order 
of say 50 comes from the n th strip itself since most absorption occurs very 
close to the point of origin to the photons. Near the origin, it is 
desirable to approximate pure absorption by a linear function rather than 
by an average value. Hence, for the n*"* 1 zone itself, we have used the 

transport function * n the integrand. In contrast to the dominance of the 

til £ 

n term in E , we were able to ignore M and hence the function as well 
n nn i 

in the evaluation of the transport integral. 

Evaluation of the integral in Equation (A27) with n = m yields the 


value 


_ ^ At At 

2[£ 2 (0,0, H n )- £ 2 <-jS, (A28) 

For the other zones it is sufficient to consider t' constant in each zone 
and use the form ^ ^ in the integrands . This leads to 



m + 1 


T 


m 
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1* 


t ) - (At , 
nm ^1 nm* 


t ) 

nm 


4X 

(A29) 


Thus 


E 
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1 - 2[,f 2 (0,0,a t n )- £ 2 (. 


At At 
n n 


n 


2 »“t J 


m 


n 


|<?. (At . t ) - f. (A t , t ) 
1 1 n,m + X’ nm 1 n,m’ nm' 


(A30) 


It is interesting to note that the addition of pure absorption actuaXXy 
increases the "escape" probabiXity since by definition, escape here refers 
to loss from the medium without a resonance absorption. 

FinaXXy the extension of this treatment to the three Xines at 

G 

X302, 1304 and I306A is simpXe. The intensity becomes 


4irl 

n 


X 

E 


W 





nm n,m + I 
- w i > 


CA31) 


while the transport solution becomes 


S on h 2> * W E n + 


kx 


m 


2 . * •'* \ 


nm 


X k 


(A32) 


where S and S are expressed in terms of depth t„. The escape function in 
o ^ 

Equation (A32) is now the sum of probabiXities of escape in each of the 
three Xines, i.e. , 


E 

n 


X- E { gj » 2[(f 2 (0,0,a t n f)-f 2 (^,^,a^)] 

. „ , \£ n,m + I nm N . . n,m 

+ z g. ICiCAt.’ ,t ) - £ ,(at ’ ,t )|} 

m 4 n ** 3 J 


(A33) 



k2 


In summary we have reduced the transport integral and the escape 

function to expressions involving the transport functions W (Equation 

nm 

(A2)) , £ (Equation (A8)), (Equation (A19)) and (Equation 

(A21)), ye considered four cases for the transport integral: 

1) No pure absorption 

a) S(t) constant over a given zone; Equation (A12) 

b) S(t) linear over a given zone; Equation (A13) 

2) Pure absorption included 

a) S(t) constant; Equation (A22) 

b) S(t) linear; Equation (A26) 

Corresponding to the transport integral are two cases for the escape 
function: 

1) No pure absorption; Equation (A9) 

2) Pure absorption included: Equation (A30). 
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FIGURE CAPTIONS 

1. Directional and frequency dependent CFR intensity profiles. Because of 
considerably greater optical depths along the line of sight in the 
downward direction (9 = 180*) the corresponding profiles at large x 
are greater than those for the upward observing direction. 

2. Directional and frequency dependent CFR intensity profiles near the 
upper boundary. 

3. Smoothed experimental 1304 intensity height profile including the low 
altitude data of Fastie et al. (1964) and high altitude data of Fastie 
(1968) . 

4. Behavior of S^) near the upper boundary for models with exospheric 
temperatures of 850* and 1000*. The dashed curves represent S^) for 
the case which ignores attenuation from the point of emission to the 
point of observation. 

5. Total source functions S(z) for the 850° exospheric temperature model. 
The solutions show the effect of changing T_. In each case T_ = 300°. 

6. Initial source functions Sq(z) corresponding to the S(z)’s of Fig. 5. 

7. Total source functions S(z) for the 850* exospheric temperature model. 
The solutions show the effect of changing T^ from 300* in model 3C to 
850® in model 8C. In each case, ^ = 1000*. 

8. Initial source functions S q (z) corresponding to the S(z)'s of Fig. 7. 

9. Initial source functions S q (z) of Fig. (6) and(8) shown for comparison. 
Included on the abscissa is the temperature profile of the 850* exo- 
spheric temperature model. 

10. Initial source functions for the 1000® exospheric temperature model and 
f - .021 shown for comparison among themselves and to compare with those 
of Fig. 9 which pertain to f = .046. 
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11, The initial source functions S ^^ (z) and g ^^(z) for the 850° and 

o o 

1000° exospheric temperature models extended to high altitudes and the 

solar resonance production rate determined explicitly by Equation (39). 

8C 

12. Initial source function S q (z) and production rates determined for 

photoelectron impact excitation on 01. All solutions refer to a solar 

zenith angle of 60°. The 850° model was used to determine S (z) as well 

o 

as Stewart's solution while the 1000° model was used to determine Prasad's 
solution. 
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Table 2. Boltzmann factors for the triplet levels of the ground state 
of 01. 


1 


200 

.812 

500 

.674 

1000 

.617 

09 

.555 


JL 

-2 

156 

.032 

256 

.070 

294 

.089 

333 

.111 
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Table 3. 850 9 Exospheric Temperature Model 


;(Kta) 

T(°K) 

n(0) 

n (0g) 

100 

208 

5.00(11) 

1.99(12) 

120 

355 

7.60(10) 

7.50(10) 

150 

643 

1.43(10) 

4.83(9) 

200 

801 

3.42(9) 

3.44(8) 

250 

838 

1.12(9) 

3.87(7) 

300 

847 

3.97(8) 

4.90(6) 

350 

849 

1.45(8) 

6.52(5) 

400 

850 

5.37(7) 

8.99(4) 

450 

850 

2.05(7) 

1.15(4) 

500 

850 

7.74(6) 

1.87(3) 

600 

850 

1.18(6) 

4.33(1) 

700 

850 

1.76(5) 

1.00(0) 


^2 

^l 

12 

7.19(4) 

1.92(4) 

4.63(3) 

1.58(4) 

6.20(3) 

1.73(3) 

5.25(3) 

2.32(3) 

6.79(2) 

1.57(3) 

7.17(2) 

2.13(2) 

5.47(2) 

2.51(2) 

7.44(1) 

2.00(2) 

9.17(1) 

2.72(1) 

7.41(1) 

3.40(1) 

1.01(1) 

2.78(1) 

1.28(1) 

3.80(0) 

1.07(1) 

4.90(0) 

1.46(0) 

4.12(0) 

1.89(0) 

5.6l(-l) 

6.24(-l) 

2.87(-l) 

8.50(-2) 

9. 77 (-2) 

4. 48 (-2) 

1*33 (-2) 
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Table 4. 1000° Exospheric Temperature Model 


(Km) 

(T°K) 

n(o) 

n (° 2 ) 

h 



100 

220 

2.83(11) 

2.90(12) 

3.26(4) 

1.04(4) 

2.58(3) 

120 

347 

9.65(10) 

1.14(11) 

1,04(4) 

4.08(3) 

1.14(3) 

150 

661 

1.64(10) 

6.52(9) 

3.35(3) 

1.52(3) 

4.50(2) 

200 

893 

3.93(9) 

5.04(8) 

1.08(3) 

5.12(2) 

1.54(2) 

250 

990 

1.41(9) 

7.10(7) 

4.32(2) 

2.06(2) 

6.24(1) 

300 

1000 

5.80(8) 

1.22(7) 

1.83(2) 

8.74(1) 

2.64(1) 

350 

1000 

2.46(8) 

2.19(6) 

7.89(1) 

3.77(1) 

1.14(1) 

400 

1000 

1.05(8) 

4.07(5) 

3,45(1) 

1.65(1) 

4.98(0) 

500 

1000 

2.03(7) 

1.52(4) 

6.84(0) 

3.27(0) 

9.876-1) 

600 

1000 

4.13(6) 

6.20(2) 

1.42(0) 

6.79(-l) 

2.05(-l) 

700 

1000 

8.72(5) 

2.78(1) 

3.04(-l) 

1.45 (-1) 

4.39(-2) 
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Table 5 « Models used to obtain transport solutions 


for 

the f values .021 and .046 


f * .021 

f = *oU6 

T ex * 1000 ° 

T ejt - 850 " 

T ex = 1000° 

Sir v° K > v° k > 

Sir V K > V’ K > 

SL T B ( " K ) 

3 V(a) 300 variable 

30(a) 300 1000 

10 V(a) 1000 variable 

10C (a) 1000 1000 

3 Vfb) 300 variable 

3C(b) 300 1000 

8v(b) 850 variable 

8c(b) 850 1000° 

10 C(c) 1000 1000 


850 variable 
850 1000° 
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